The  Ohio  State  University 


TRANSIENT  CURRENT  DENSITY  WAYEFORMS  ON 
A  PERFECTLY  CONDUCTING  SPHERE 

E.M.  Kennaugh 
D.L.  Moffatt 


The  Ohio  State  University 

ElectroScience  Laboratory 

Department  of  Electrical  Enginer-ing 
Columbus,  Ohio  43212 

Technical  Report  710816-13 
Contract  N00014-78-C-0049 
October  1983 


o 

o 

.  ■ 

Department  of  the  Navy 

Office  of  Naval  Research 

800  N.  Quincy  Street 

d 

Arlington,  Virginia  22217 

DISTRIBUTION  STATEMENT  A 

Approved  fox  public  release! 
Distribution  Unlimited 


s 


DTIC 

EL.ECTE 
DEC  1 5  1983 

D 


83  It  08  08+ 


■  ■  as*  •  bv  , 


NOTICES 


When  Government  drawings,  specifications,  or  other  data  are 
used  for  any  purpose  other  than  In  connection  with  a  definitely 
related  Government  procurement  operation,  the  United  States 
Government  thereby  Incurs  no  responsibility  nor  any  obligation 
whatsoever,  and  the  fact  that  the  Government  may  have  formulated, 
furnished,  or  In  any  way  supplied  the  said  drawings,  specifications, 
or  other  data,  is  not  to  be  regarded  by  Implication  or  otherwise  as 
In  any  manner  licensing  the  holder  or  any  other  person  or  corporation, 
or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  Invention  that  may  In  any  way  be  related  thereto. 


I 


SOI  77  -Itll 

KNOUT  DOCUMENTATION  j  i.  MW  ■«.  |* 

_ «si _  /vyA/zc  zo>?~ _ i _ 

4.  TWO  and  KuMiTI* 

TRANSIENT  CURRENT  DENSITY  WAVEFORMS  ON  A 
PERFECTLY  CONDUCTING  SPHERE 

*  Kennaugh  and  D.L.  Moffatt 

t.  Parfcrmint  Of*»mi»ti«n  Nam*  and  Addfata 

The  Ohio  State  University  ElectroScience  Laboratory 
Department  of  Electrical  Engineering 
Columbus,  Ohio  43212 

It.  Bpontaring  Organisation  Nam*  and  Addfatt 

Department  of  the  Navy,  Office  of  Naval  Research 
800  North  Quincy  Street  (614  A:CEA) 

Arlington,  Virginia  22217 

If.  Suppramantary  Nout 


I.  baeip.*nt  a  Accaoion  No 

jfT  Rape*  Dal* 

January  1983 


B.  Partormin*  Organnation  Rapt.  No 

F.SL  710816-13 

10.  Rrotoet/Taak/Werli  Unit  No. 


)1.  Contract(C)  at  Crant(G)  No 

to 

,C)  N00014-78-C-0049 

If.  Typo  ot  Raport  4  Ranatl  Cavarad 

Technical  Report 


—  —  .  -  -  — — 

Abstract  (Limit:  700  words) 

^The  time-dependent  surface  current  density  waveforms  at  various  points  on  the 
surface  of  a  perfectly  conducting  sphere  Illuminated  by  a  plane  electromagnetic  wave  ar< 
presented.  The  incident  plane  wave  has  an  impulsive  or  shock-type  time  dependence. 
Currept  density  waveforms  In  the  E-plane  (strong  creeping  wave)  and  H-plane  on  both  the 
t""11t*-and  shadowed  sides  of  the  sphere  are  shown.  On  the  illuminated  side  of  the 
sphere,  removal  of  the  Klrchhoff  or  physical  optics  approximation  for  the  surface 
current  density  permits  a  detailed  examination  of  the  deficiencies  of  this  approximate 
for  short  times  (high  frequencies).  From  these  results  the  form  of  a  first  order 
correction  to  the  Klrchhoff  current  Is  suggested.  On  the  shadowed  side  of  the  sphere, 
the  changing  character  of  the  current  waveform  is  noted. 
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a  perfectly  conducting  sphere.  Illuminated  side  of  sphere. 

2.  E-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 

3.  H-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 

4.  E-plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  current  density  has  been  removed. 

Illuminated  side  of  sphere. 

5.  E-plane  (solid)  and  H-plane  (dashed)  error  for  specular 
term  (exact  minus  physical  optics)  for  surface  current 
density  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere. 

6.  E-plane  (solid)  and  H-plane  (dashed)  error  for  specular 
term  (exact  minus  corrected  physical  optics,  value  of 

K  constant  Is  0.2)  for  surface  current  density  induced 
by  an  Impulsive  plane  wave  on  a  perfectly  conducting 
sphere. 

7.  E-plane  (solid)  and  H-plane  (dashed)  error  for  specular 
term  (exact  minus  corrected  physical  optics,  value  of 

K  constant  Is  0.1)  for  surface  current  density  Induced 
by  an  Impulsive  plane  wave  on  a  perfectly  conducting 
sphere. 
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I.  INTRODUCTION 


Development  of  the  Impulse  response  concept  In  3-dlmenslonal 
electromagnetic  scattering  problems  [1,2]  was  materially  aided  by  a 
Fourier  synthesis  procedure  whereby  harmonically  related  complex 
scattering  data  were  used  to  synthesize  approximate  far-zone  canonical 
response  waveforms  [2, 3, 4, 5].  The  same  procedure  permits  one  to  generate 
the  time-dependent  surface  current  density  waveforms  on  a  scatterer  or 
radiator.  A  comparison  of  Fourier  synthesis  and  space-time  Integral 
equation  calculations  has  been  made[6],  Both  procedures  have  an  Inherent 
resolution  limitation  dictated  either  by  &  finite  summation  (Fourier 
synthesis)  or  a  Gaussian  pulse  excitation  (space-time  Integral  equation), 
but  it  Is  not  correct  to  suggest  that  an  Infinite  summation  [7]  or 
equivalently  an  Impulsive  excitation  must  be  used.  However,  some  care 
must  be  exercised  In  Interpreting  singularities  and/or  jump 
discontinuities  particularly  If  they  occur  simultaneously. 

Impetus  for  this  report  stems  from  a  depolarization  correction  to 
physical  optics  derived  by  Bennett  [8]  and  utilized  by  Boerner  [9].  We 
show  that  even  for  an  object  where  no  depolarization  occurs  (the 
conducting  sphere),  the  form  of  a  basic  correction  to  the  Klrchhoff 
current  can  be  deduced.  It  Is  suggested  that  similar  first  order 
corrections  will  Improve  the  Klrchhoff  approximation  for  a  number  of 
different  object  geometries. 


A.  Surface  Current  Density  Waveform 


Hie  series  calculations  of  the  tangential  magnetic  field  and  radial 
electric  field  at  the  surface  of  a  conducting  sphere  of  radius  a  Immersed 
In  an  Incident  plane  wave  have  been  made  for  sphere  cl rcumf erences  In 
wavelengths  of  0.2  (0.2)20.0.  Locations  on  the  sphere  surface  In  both 
the  E-plane  and  H-plane  for  theta  angles  of  0  (15)  180  degrees  were  used 
with  0  degrees  crr^spondlng  to  the  specular  point.  The  real, 
time-dependent  turf ace  current  density  waveforms  synthesized  from  these 
data  are  shown  In  Figures  1,  2  and  3.  Figure  1  combines  the  E-plane 
(positive)  and  H-plane  (negative)  waveforms  on  the  Illuminated  side  of 
the  sphere.  The  abscissa  scale  Is  In  units  of  transit  time  for  the 
sphere  diameter.  On  the  Illuminated  side  of  the  sphere,  the  weighted 
quasi -Impulses  at  the  origin  of  each  waveform  Illustrate  the  resolution 
obtained  (approximately  o.75a).  On  the  Illuminated  side  of  the  sphere, 
the  H-plane  waveforms  have  a  sign  reversal.  On  the  Illuminated  side  of 
the  sphere,  the  current  density  Is  primarily  Impulsive;  details  for  short 
times  are  best  seen  after  the  Klrchhoff  approximation  for  the  current  Is 
removed.  There  Is  clear  evidence  of  the  creeping  wave  contributions. 

The  surface  current  density  waveforms  on  the  shadowed  side  of  the  sphere 
are  shown  In  Figure  2  (E-plane)  and  Figure  3  (H-plane).  Note  that  the 
ordinate  scales  for  H-plane  are  not  constant.  On  the  shadowed  side  of 
the  sphere  note  that  the  character  of  the  current  density  at  the  onset  of 
the  waveform  Is  still  apparently  Impulsive  near  the  shadow  boundary 
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(theta  less  than  120  degrees)  In  the  H-plane  but  not  In  the  E-plane. 

There  Is  a  strong  creeping  wave  contribution  In  the  E-plane  and  some 
evidence  of  a  very  weak  creeping  wave  contribution  In  the  H-plane. 

In  Figure  4  the  surface  current  density  waveforms  on  the  Illuminated 
side  of  the  sphere  are  shown  after  the  Klrchhoff  approximation  to  the 
current  density  has  been  removed.  The  H-plane  waveforms  are  dashed.  It 
Is  very  evident  that  as  the  observation  point  progresses  from  the 
specular  point  to  the  shadow  boundary,  the  Klrchhoff  approximation  becomes 
progressively  too  large  In  the  E-plane  and  too  small  In  the  H-plane.  At 
a  given  angle,  the  magnitude  of  the  errors  In  the  E-plane  and  H-plane  are 
approximately  equal.  This  would  seem  to  Imply  that  on  a  cut  at  45 
degrees  from  the  E-plane  and  H-plane  the  Klrchhoff  approximation  to  the 
current  may  be  quite  good  even  near  the  shadow  boundary.  In  the  E-plane 
at  theta  equal  60  degrees  and  75  degrees  one  can  also  clearly  discern  two 
distinct  creeping  wave  contributions,  both  launched  from  the  shadow 
boundary  and  delayed  from  the  waveform  onset  by  different  travel  times  on 
the  sphere  surface. 

B.  Correctl on 

Based  on  the  results  shown  for  the  conducting  sphere,  the 
following  formulas  for  correcting  the  physical  optics  estimates  for  the 
Induced  currents  are  suggested  for  smooth  (no  edge)  scatterers.  In  the 
E-plane  modify  the  physical  optics  currents  by  the  multiplicative  factor 
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[1  -  K  sins] 


(1) 


In  the  H-plane  modify  the  physical  optics  currents  by  the  multiplicative 
factor 


[1  ♦  K  ±12®_]  ,  (2) 

cos  6 

where  e  Is  the  angle  between  the  outward  normal  to  the  surface  of 
scatterer  and  the  backscatter  direction.  The  H-plane  correction  becomes 
Infinite  at  the  shadow  boundary  (8«90*)  but  when  multiplied  by  the 
physical  optics  current  estimate  (cos 8  factor)  yields  a  finite  current 
estimate.  The  specular  error  (exact  minus  physical  optics)  before  cor¬ 
rection  Is  shown  In  Figure  5.  Figures  6  and  7  show  the  effect  of  the 
correction  term  for  values  of  the  constant  K  of  0.2  and  0.1,  respective¬ 
ly.  For  the  simple  form  of  the  correction  suggested,  a  value  K  near 
0.1  appears  optimum.  Other  forms  for  the  correction,  possibly  exponen¬ 
tial,  may  suggest  themselves  to  the  reader.  Our  purpose  here,  however, 
was  to  show  that  simple  corrections  are  possible. 

C.  Conclusions 

It  has  been  demonstrated  that  relatively  simple  formulas  can  be 
developed  for  correcting  the  p^  si cal  optics  estimates  of  the  current 
densities  Induced  on  a  conducting  spherical  scatterer.  It  has  also  been 
suggested  that  similar  simple  corrections  could  be  made  for  o,her 
"smooth"  scatterers.  It  Is  noted  that  the  corrections  Indicated  here  are 
not  related  to  depolarization  properties  of  the  scatterer. 


The  surface  current  density  waveforms  on  the  shadowed  side  of  the 
sphere  have  been  Included  here  for  the  possible  benefit  of  other 
researchers.  The  physical  optics  estimate  has  long  been  a  favorite 
direct  scattering  solution  for  Inverse  scattering  because  the  estimate 
can  be  generally  related  to  a  target  cross  sectional  area  function  on  the 
Illuminated  side  of  the  scatterer.  A  similar  relationship  (general)  for 
the  shadowed  side  of  the  scatterer  would  greatly  Improve  Imaging 
techniques  based  on  the  physical  optics  estimate. 


a.  e  *  0#  (specular) 


Figure  1.  E-plane  (positive)  and  H-plane  (negative)  surface  current 
density  waveforms  Induced  by  an  Impulsive  plane  wave  on 
a  perfectly  conducting  sphere.  Illuminated  side  of  sphere. 
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Figure  1.  E-plane  (positive)  and  H-plane  (negative)  surface  current 
density  waveforms  Induced  by  an  Impulsive  plane  wave  on 
a  perfectly  conducting  sphere.  Illuminated  side  of  sphere. 
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Figure  1. 


E-plane  (positive)  and  H-plane  (negative)  surface  current 
density  waveforms  Induced  by  an  Impulsive  plane  wave  on 
a  perfectly  conducting  sphere.  Illuminated  side  of  sphere 
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Figure  1.  E-plane  (positive)  and  H-plane  (negative)  surface  current 
density  waveforms  Induced  by  an  Impulsive  plane  wave  on  a 
perfectly  conducting  sphere.  Illuminated  side  of  sphere. 
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Figure  1.  E-plane  (positive)  and  H-plane  (negative)  surface  current 
density  waveforms  Induced  by  an  Impulsive  plane  wave  on  a 
perfectly  conducting  sphere.  Illuminated  side  of  sphere. 
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Figure  1.  E-plane  (positive)  and  H-plane  (negative)  surface  current 
density  waveforms  Induced  by  an  Impulsive  plane  wave  on  a 
perfectly  conducting  sphere.  Illuminated  side  of  sphere. 


11 


* 


g.  e  ■  90* 


Figure  1.  E-plane  (positive)  and  H-plane  (negative)  surface  current 
density  waveforms  Induced  by  an  Impulsive  plane  wave  on  a 
perfectly  conducting  sphere.  Illuminated  side  of  sphere. 
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Figure  t.  E-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  vave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 


E-plane  surface  current  density  waveforms  Induced  by  an 
impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sohere. 


c.  e  «  135* 


Figure  2.  E-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 


Figure  2.  E -plane  surface  current  density  waveforms  Induced  by  an 

Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 


e.  8  »  165* 


Figure  2.  E-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 
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f.  8  -  180# 


Figure  2.  E-plane  surface  current  dens It/  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 
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Figure  3.  H-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 
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b.  e  ■  120 


Figure  3. 


H-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 
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c.  e  -  135* 


Figure  3.  H-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 


Figure  3.  H-plane  surface  current  density  waveforms  Induced  by  an 
Impulsive  plane  wave  on  a  perfectly  conducting  sphere. 
Shadowed  side  of  sphere. 
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Figure  4.  E-plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  current  density  has  been  removed. 

Illuminated  side  of  sphere. 
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Figure  4. 


b.  e  -  15* 


E -plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  currrent  density  has  been  removed. 
Illuminated  side  of  sphere. 
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c.  8  ■  30# 


Fig-.  -a  4.  E-plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  current  density  has  been  removed. 

Illuminated  side  of  sphere. 
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Figure  4.  E-plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  current  density  has  been  removed. 

Illuminated  side  of  sphere. 


Figure  4.  E-plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  current  density  has  been  removed. 

Illuminated  side  of  sphere. 
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Figure  4. 
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E-plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  current  density  has  been  removed. 

Illuminated  side  of  sphere. 
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g.  e  «  90* 


Figure  4.  E-plane  (solid)  and  H-plane  (dashed)  surface  current  density 
waveforms  Induced  by  an  Impulsive  plane  wave  on  a  perfectly 
conducting  sphere  after  the  physical  optics  approximation 
for  the  surface  current  density  has  been  removed. 

Illuminated  side  of  sphere. 


- E- PLANE 


Figure  5.  E-plane  (solid)  and  H-plane  (dashed)  error  for  specular 
ten*  (exact  «rlnus  physical  optics)  for  surface  current 
density  Induced  by  an  Inpulslve  plane  wave  on  a  perfectly 
conducting  sphere. 
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- E- PLANE 

- H- PLANE 


Figure  6.  E-plane  (solid)  and  H-plane  (dashed)  error  for  specular 
term  (exact  minus  corrected  physical  optics,  value  of 
K  constant  Is  0.2)  for  surface  current  density  Induced 
by  an  Impulsive  plane  wave  on  a  perfectly  conducting 
sphere. 
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- E- PLANE 

-  H-PLANE 


Figure  7.  E-plane  (solid)  and  H-plane  (dashed)  error  for  specular 
term  (exact  minus  corrected  physical  optics,  value  of 
K  constant  Is  0.2)  for  surface  current  density  Induced 
by  an  Impulsive  plane  wave  on  a  perfectly  conducting 
sphere. 
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